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ABSTRACT 

We show that atomic alignment presents a reliable way to study topology of astrophysical magnetic 
fields. The effect of atomic alignment arises from modulation of the relative population of the sublevels 
... of atomic ground state pumped by anisotropic radiation flux. As such aligned atoms precess in the 

' external magnetic field and this affects the properties of the polarized radiation arising from both 

I scattering and absorption by the atoms. As the result the polarizations of emission and absorption 

. lines depend on the 3D geometry of the magnetic field as well as the direction and anisotropy of 

incident radiation. We consider a subset of astrophysically important atoms with hyperfine structure. 
O \ For emission lines we obtain the dependencies of the direction of linear polarization on the directions 

O , of magnetic field and the incident pumping radiation. For absorption lines we establish when the 

polarization is perpendicular and parallel to magnetic field. For both emission and absorption lines 
we find the dependence on the degree of polarization on the 3D geometry of magnetic field. We claim 
that atomic alignment provides a unique tool to study magnetic fields in circumstellar regions, AGN, 
interplanetary and interstellar medium. This tool allows studying of 3D topology of magnetic fields 
\ and establish other important astrophysical parameters. We consider polarization arising from both 

^ . atoms in the steady state and also as they undergo individual scattering of photons. We exemplify 

' the utility of atomic alignment for studies of astrophysical magnetic fields by considering a case of Na 

' alignment in a comet wake. 

Subject headings: ISM: atomic processes — magnetic fields — polarization 



1. INTRODUCTION 



O ^ 

(-H \ Magnetic fields play essential roles in many astrophysical circumstances. Incompatible with its importance, our 
O • knowledge about astrophysical magnetic field is very much limited. Exploring new tools to measure magnetic field is 
I ' thus extremely important. 

O In our previous paper (Yan & Lazarian 2006, henceforth Paper I), we discussed how the alignment of fine structure 
^ atoms and ions can be used to detect 3D orientation of magnetic field in diffuse medium. As in Paper I, for the sake 
^ ■ of simplicity, we shall term the alignment of both atoms and ions atomic alignment. In this paper, we shall discuss 
. . \ atomic alignment within hyperfine structures. In fact, historically optical pumping of atoms with hyperfine structure 
K*" ■ atoms, e.g., alkali atoms. He, Hg at al. (Happer 1972) were first studied in laboratory in relation with early day 
\^ maser research. This effect was noticed and made use of for the interstellar case by Varshalovich (1968). And in a 
rS \ subsequent paper by Varshalovich (1971), it was first pointed out that the dependence of atomic alignment on the 
■ direction of magnetic field can be used to detect magnetic field in space. However, the study did not provide either 
■ ■ ■ ' detailed treatment of the effect or quantitative predictions. 

Atomic alignment has been addressed also by solar researchers. The research into emission line polarimetry resulted 
in important change of the views on solar chromosphere (see Landi Degl'Innocenti 1983, 1984, 1998, Stcnflo & Keller 
1997, Trujillo Bueno & Landi Degl'Innocenti 1997, Trujillo Bueno 1999, Trujillo Bueno et al. 2002, Manso Sainz & 
Trujillo Bueno 2003). However, they dealt with emission of atoms in a very different setting. Similar to Paper I 
below we concentrate on the weak field regime, in which it is the atoms at ground level that are repopulated due to 
magnetic precession, while the Hanle effect that the aforementioned works deal with is negligible. The closest to our 
study is regime discussed in the work by Landolfi & Landi Degl'Innocenti (1986), who considered an idealized two-level 
fine structure atom in a very restricted geometry of observations, namely, when the magnetic field is along of line of 
sight and perpendicular to the incident light. As we discussed in Paper I, that dealt with fine structure atoms, these 
restrictions did not allow to use this study to predict the directions of astrophysical magnetic fields from polarimetric 
observations. 

In this paper we consider atomic alignment for atoms with hyperfine structure and provide quantitative predictions for 
both absorption and emission lines. To exemplify the processes of alignment we perform calculations for astrophysically 
important species. 

Studies of magnetic field topology using atomic alignment are complementary to the studies of magnetic fields 
using aligned dust. Similar to the case of interstellar dust, the rapid precession of atoms in magnetic field makes 
the direction of polarization sensitive to the direction of underlying magnetic field (see Lazarian 2003 for a review). 
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Fig. 1. — Left: a toy model to illustrate how atoms are aligned by anisotropic light. Atoms accumulate in the ground sublevel M = as 
radiation removes atoms from the ground states M = 1 and M = —1; Right: Radiation geometry and the polarization vectors in a given 
coordinate system. O is the direction of radiation, x the positional angle of a linear polarization. 

However, as the precession of magnetic moments of atoms is much faster than the precession of magnetic moments of 
grains, atoms can reflect much more rapid variations of magnetic field. More importantly, alignable atoms and ions 
can reflect magnetic fields in the environments where either the properties of dust change or the dust cannot survive. 
This opens wide avenues for the magnetic field research in circumstellar regions interstellar medium, interplanetary 
medium, intracluster medium, AGN, etc. In addition, the polarization caused by atomic alignment is sensitive to the 
3D direction of magnetic field. This information is not available by any other technique that are available for studies 
of magnetic field in diffuse gas. 

In what follows we formulate the conditions for atomic alignment of hyperfine species in §2 and present our formalism 
for treating of atomic alignment and optical pumping in §3. In §4 we use Nal and KI to discuss the details of practical 
calculations of polarization of emission lines arising from atomic alignment of species with hyperfine splitting. A 
nonequilibrium case is considered for Na I in §5 and it is shown how alignment and polarization changes with the 
number of scattering events. The results are then applied to turbulent comet wake in §6. We show how the change 
of magnetic field direction infiuences the polarization of scattered Sodium D lines from the wake which can be used 
to ground-based studies of interplanetary turbulence. In §7, we consider the alignment of neutral hydrogen, N V and 
P V and the resulting polarizations of Lyman a. In addition, we also discuss their implications for HI 21cm and N V 
hyperfine radio lines. More complicated atomic species are considered in §8, where we show that absorption from N 
I atoms is polarized even for unresolved hyperfine multiplet. We show how hyperfine structure changes its alignment 
compared to S II, which has the same electron configuration, but without nuclear spin. In §9, we discuss how the 
average along the line of sight affects the results. The discussion and the summary are provided in, respectively, §10 
and §11. 



We discussed atomic alignment in Paper I. Atomic alignment is caused by the anisotropic deposition of angular 
momentum from photons. As illustrated by the toy model shown in FigOe/i, absorption from Alp = in the ground 
level is impossible owing to conservation of angular momentum. The differential absorptions for a realistic example 
of Nal will be provided later in §4.1. As the result, atoms scattering the radiation from a light beam are aligned in 
terms of their angular momentum. To have the alignment of ground state, the atom should have non-zero (> 1) total 
angular momentum on its ground state to enable various projection of atomic angular momentum. 

For atoms with nuclear spin, hyperfine structure need to be taken into account. It is the total angular momentum 
J + 1 = F, the vector summation of electron angular momentum J and nuclear spin I, that should be considered. 
Alkali atoms thus are alignable with nuclear spins added. It is true that direct interaction between the nucleus and 
the radiation field is negligible. However, for resonant lines, the hyperfine interactions cause substantial precession 
of electron angular momentum J about the total angular momentum F before spontaneous emission. Therefore total 
angular momentum should be considered and the FMp basis must be adopted (Walkup, Migdall & Pritchard 1982). 
For alkali-like atoms, hyperfine structure should be invoked to allow enough degrees of freedom to harbor the alignment 
and to induce the corresponding polarization. 

In order for atoms to be aligned, the collisional rate should not be too high. In fact, as disalignment of the ground 
state requires spin exchange (or fiips), it is less efficient than one can naively imagine. The calculations by Hawkins 
(1955) show that to disalign sodium one requires more than 10 collisions with electrons and experimental data by 
Kastler (1956) support this. This reduced sensitivity of aligned atoms to disorienting collisions makes the effect 
important for various astrophysical environments. 

Magnetic field would mix up different M states. However, it is clear that the randomization in this situation is not 



2. CONDITIONS FOR ATOMIC ALIGNMENT IN THE PRESENCE OF HYPERFINE STRUCTURE 
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TABLE 1 

Selected alignable atomic species and corresponding transitions. Note only lines above 912A are listed. Hyperfine 

STRUCTURE SHOULD BE CONSIDERED FOR THESE ELEMENTS WITH NUCLEAR SPINS. NOTE THAT ALKALI-LIKE SPECIES ARE NOT ALIGNABLE 

WITHOUT HY PERFINE STRUCTURE TAKEN INTO ACCOUNT (SEE TEXT). DATA ARE TAKEN FROM THE ATOMIC LiNE LiST 

|hTTP: / /www. PA. UKY.EDU/~PET ER/ATOMIO/ The last column gives the maximum polarizations FOR EMISSION BY ALKALI SPECIES 
AND ABSORPTION BY T^l; FOR OTHER SPECIES, "CNA" STANDS FOR CURRENTLY NOT AVAILABLE. 

complete and the residual alignment and resulting polarizations reflects the magnetic field direction. Magnetic mixing 
happens if the angular momentum precession rate is higher than the rate of the excitation of atoms from the ground 
state, which is true for many astrophysical conditions. 

3. PHYSICS OF ATOMIC ALIGNMENT 

We consider a realistic atomic system, which can have multiple upper levels J„ , and lower levels Fi . When such 
an atomic system interacts with resonant radiation, there will be a photoexcitation foUowe^l by spontaneous emission. 
We describe the atomic occupation and radiation field by irreducible density matrices p^, Jq (see ApplBj. Unlike fine 
structure levels, the hyperfine separation vp^^pi^ = \E(Fu) — E{F^)]/h is comparable to the natural line-width A, and 
therefore the coherence between hyperfine levels must be taken into account on the upper level. The density matrix 
of the upper state is thus There is no coherence on the ground state, and its occupation thus can be 

characterized by p^{Fi). These density matrices are determined by the balance of the three processes: absorption (at 
a rate ~ B/„ Jg), emission (at a rate ~A) and magnetic precession (at a rate ~ vl) among the sublevels of the state. 
The statistical equilibrium equations for hyperfine transitions can be extrapolated from the case of fine transitions 
(Paper I, see also Landi Degl'Innocenti & Landolfi 2004). The formalism for hyperfine transitions can be obtained by 
replacing J , Mj with F , Mp and taking into account those additional factors in Eq. (|A2D . 

p^q{Fu,Fl^) + 2niVLguqp':,{Fu,Fi) + 2niVF^F,/^{Fu, F^) = -A( J„ ^ Ji)pl{Fu,Fi) 

+ [Ji] J2 (^kk'Pk'Bi^J^ + rkk'Bi.J^)p'^^,iFl) (1) 



p^giFi) + 2TTi,,Lgiqp'qiFi)^ Pk[Ju]A{Ju^ Ji)pl{F^,Fi)- ihk'BiJ° + Skk'BiJ^)p'l^,iFi), (2) 



where 



« ^ {pJll} K.Ki'" {i?;; }} , (3) 
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m.' = {-lf+'''i3[k,k',2]y/'[F^,F:,]'/^[F[ 




Ju Jl 1 

F! FL I 



(4) 



J/ J/ 2 
1 1 Ju 



JiJi2 
FiFiI 



(5) 



In Eqs(I3][5]), the matrices with big "()" are 3-j symbol and the matrices with big "{ }" represents the 6-j or 9-j symbol, 
depending on the size of the matrix. They are also known as Wigner coefficients (see Cowan 1981, Paper I). Throughout 
this paper, we define [j] = 2j + 1, which means [F, F'] = {2F + 1)(2F' + 1). The evolution of upper state iPq{Fu,Fl^)) 
is represented by Eq.([T]) and the ground state {Pg{Fi)) is described by Eq.Q. The second terms on the left side of 
Eq. (|ll2p represent mixing by magnetic field, where and gi are the Lande factors for the upper and ground level. For 
the upper level, the mixing ~ VLguPqiFu, Fu) is much slower than the emission ^ Ap'^{Fu, F^^) and is thus neghgible 
as we consider a regime where the magnetic field is much weaker than Hanle field^. {i^l ^ -A/gu)- The third term on 
the left side of Eq.([T|) gives a measure of coherence of two hyperfine levels. It is easy to see if Vp^ p^ ^ A, the Einstein 

emission coefficient, the coherence component of the density matrix p'^{Fu, F^J) would be zeros. The two terms on 
the right side of Eq. ([Tl ^ are due to spontaneous emissions and the excitations from ground level. Transitions to all 
upper states are taken into account by summing over J„ , F^ in Eq.Q. Vice versa, for an upper level, transitions to 
all ground sublevels {Fi) are summed up in Eq.([T|). The excitation is proportional to 



J, 



K 



4^^ 

1=0 



(6) 



which is the radiation tensor of the incoming light averaged over the whole solid angle and line profile ^(i^ — vq) 
Si = [I, Q, U, V] represent Stokes parameters. The unit radiation tensors jQ{i,Vl) are given by: 

(1 - 1.5 sin^ e) 



J'|2(^,^!) = V3e±2*■^ 



For unpolarized point source from 



-3/2sin^6' 







sin^ e/A 
-(1 + CC 



^0)/4 



=F«cos6'/2 
ir), the radiation tensor is then: 



=Fsin26'/4\ 
T sin 26*74 
-zsin6'/2 / 



JO J2 = J^(2-3sin2 9r)F,Jl2 = sin2 0,e±*20.j j2 



sin26'^e±*'^-/* 

AW 



(7) 



(8) 



where W is the dilution factor of the radiation field, which can be divided into anisotropic part Wa and isotropic part 
Wi (Bommier & Sahal-Brechot 1978), /* is the solid-angle averaged intensity. In the case of a point source, Wi = 0. If 
Wi 0, the degree of alignment and polarization will be reduced. The solid-angle averaged intensity for a black-body 
radiation source is ^ 

Since we are interested in the regime where for the ground level the magnetic mixing is much faster than the optical 
pumping vl ^ = BJq for the ground state, magnetic coherence does not exist either. Thus there are only 
components Pq{Fi) for the ground level. Taking into account this simplification, we obtain steady state solutions by 
setting the first terms of Eq. ()l|2|l on the left side to zeros, 

NT^ 1 



l + 2TrivF^p'/A 

r u r ^ ^ 



Pk[Ju, Jl] } (Skk'Pk'BiuJo + rkk'BiuJo)Po (Fl) 



F.'h 



J^k' 



{Skk'BiuJo + Skk' BiuJ^)p'o {Fi)^0 



(10) 



Bii 
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\J 



■ F' 



^{Skk'Pk'Ja 



- Tkk 



JDpUfi) 



(11) 



F,k' 



^ For the Hanle effect to be dominant, magnetic splitting ought to be comparable to the energy width of the excited level. 
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It can be proved from Ea. (fTO| that Pq cx Jg . Eq. (fTO|) represents a set of linear equations. Considering the equation 
with fc = 0, it only includes pg' due to the triangular rule of the 3j symbol in the coefficient rkk' ■ For pg, the coefficient 
is cx Jq as the coefficient cx Jg is zero. Therefore pg cx — pg cx Jg. As a result, the dipole component of density matrix 
changes its sign at Van Vleck angle as we shall show later. This is a generic feature of atomic alignment independent 
of their specific structures of atomic levels. The corresponding emission coefficient can be extrapolated from the one 
for fine structure atoms (see Landi Degl'Innocenti 1982) by replacing (L,S,J,M) with {J, I, F, Mp), 

^ KQF^F^Fi ^ '■ " 

P^iF^,Fi)J^{z,n), (12) 

where n is the total number density of the atoms. This is the expression if we can resolve the hyperfine components 
Fi = 1 and Fi = 2. For the Dl lines of alkali atomic species, polarization will be zero otherwise. The corresponding 
emissivities in unresolved case can be found in Landolfi and Landi Degl'Innocenti (1985). Since the separations among 
the hyperfine levels on the upper state is much smaller, the absorption coefficients in this case can be obtained by 
making the analogy with the emissivities of the unresolved case, 

,.(.,17) = ^B<(.-.o)M E [i^/]^/3(-l)l-^"+^+^' P'^'f I 

KQFi III J 

]^]^|'Jp^(F,)jQ^(z,f)), (13) 
For optically thin case, the linear polarization degree and the positional angle 

1 1 



p = + uyi = ^el + ej/eo, X - - tan-i(C//Q) = - tun'^e^/e,) (14) 

(see Fig[T]); the polarization produced by absorption through optical depth t = rjod is 

Q -mdJo Vi 



It (1 - riod)IoT]Qd rjo 



U^O. (15) 



The 6-j symbol in Ea. (fT3| \ ^ r = for K — 2 and J/ < 1. This suggests that absorption is unpolarized 

y i'l Ii J J 

for atoms with J; < 1. Alkali atoms can only produce polarized emissions therefore. 

4. ALIGNMENT OF NA I, K I 

4.1. Dl and D2 lines of Na I 

The geometry of the radiation system is illustrated by Figl^l The origin of this frame is defined as the location of 
the atomic cloud. The line of sight defines z axis, and together with direction of radiation, they specify x-z plane. The 
x-y plane is thus the plane of sky. In this frame, the incident radiation is coming from (0g, 0), and the magnetic field 
is in the direction {Ob, 4>b)- 

The magnetic field is chosen as the quantization axis (z") for the atoms. Alignment shall be treated in the frame 
(see ApplDlhow these two frames are related). In this "theoretical" frame, the line of sight is in {9, tt) direction 
(i.e., the x"-z" plane is defined by the magnetic field and the line of sight, see Yig^ighi) , and the radiation source is 
directed along {9r,(t>r)- 

The ground state of Na is 2^5i/2 and the first excited states 2'^Pi/2, 2^P3/2 correspond to Dl and D2 lines respectively. 
The nuclear spin of Na is / 3/2, its total angular momentum thus can bei^ = /±J = l,2 (see Fig^ejf). According 
to the selection rule of the 3j symbol in Eq. (|Bl|) . the irreducible density tensor of the ground state Pg{Fi — 1) has 
components with fc = 1 - 1, 1, 1 + 1 = 0, 1, 2, p'^{Fi = 2) has components with fc = 2 - 2,2 - 1, ...,2 + 2 = 0, 1,2,3,4. 
For unpolarized pumping light, we only need to consider even components. For the upper level p^(-F'„), according 
to Eq.(fT2|. only the components with k — K < 2 need to be counted. For Na D2 line, the hyperfine splittings of 
upper level 2^P^^2 are comparable to their natural width: wio = 2.67, W21 = 5.27, W32 = 7.77. Thus the interference 
between levels (measured by the factor 1/(1 + 2TTii'p^^p^/A) in Ea. (|10llip must be taken into account on the upper 
level. On ground level, there is no magnetic coherence term, namely, pj^g = 0. Owing to the triangle rule of "3j" 
symbols (see Eqs llOI4(5|) . only Jq^q components appear and they are determined by the polar angle Or of the radiation 

(Eq[51 Fig[21). As a result, pg'"* are real quantities and they are independent of the azimuthal angle 0^ (see EqfTU)). 
Physically this results from fast procession around magnetic field. Their dependence on polar angle is shown in FiglHl 
To calculate the alignment for the multiplet, all transitions should be counted according to their probabilities even if 
one is interested in only one particular line. This is because they all affect the ground populations and therefore the 
degree of alignment and polarization. This means that in practical calculations, a summation should be taken over all 
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Fig. 2. — Le/t: typical astrophysical environment where atomic alignment can happen. A pumping source deposits angular momentum to 
atoms in the direction of radiation and causes differential occupations on their ground states. In a magnetized medium where the Larmor 
precession rate vj^ is larger than the photon arrival rate t^^, however, atoms are realigned with respect to magnetic field. Atomic alignment 
is then determined by the angle between the magnetic field and the pumping source. The polarization of scattered line also depends on 
the direction of line of sight, d and (or '/'n defined afterward); Middle: geometry of the observational frame. In this frame, the lino of 
sight is z axis, together with the incident light, they specify the x-z plane. Magnetic field is in {0b,<I>b) direction; Right: transformation to 
the "theoretical frame" where magnetic field defines z" axis. This can be done by two successive rotations specified by Euler angles (03, 9r) 
(see App[D]for details). The first rotation is from xyz coordinate system to x'y'z' coordinate system by an angle tpB about the z-axis, the 
second is from x'y'z' coordinate system to x"y"z" coordinate system by an angle 6 about the y'-axis. Atomic alignment and transitions are 
treated in the "theoretical" frame where the line of sight is in {9, n) direction and the incident radiation is in {9r, 4>r) direction. 

the hyperfine sublevels F„ of both the upper levels J„ — 1/2,3/2 in Eg. pTI)) . The key is the coefficients Pk,fkk'iSkk' 
which are determined by the hyperfine structure of an atomic species. By inserting the values of J„, F^, Ji, Fi, k, fc', 
K and Jq (Eq^ into Eqs.Q-dS]), we get the coefficients as given by Table 

We see from Table HTTl the coefficient Skk' =0 and this is actually true for all the alkali species. Since Skk' represents 
the differential excitation from the ground level (see Eql2]), this means alkali species are not aligned by the same 
mechanism^ (so called depopulation pumping) as illustrated by the toy model (§2). Instead they are aligned through 
another mechanism, repopulation pumping. Atoms are repopulated as a result of spontaneous decay from a polarized 
upper level (see Happer 1972). Upper level becomes polarized because of differential absorption rates to the levels 
(given by rkk', see EqU]). For instance, — —0.1042 for the absorption from i^; = 1 to F^ — 1. As the result, 
the density component of the upper level F^ = 1, Pq(u) < 0, indicating that atoms are accumulated in the sublevel 
FiMp = 1, according to the definition of irreducible tensor = [p{l, 1) — 2p(l, 0) + yo(l, —1)] (see App lBll) . 

With the coefficients Pkjfkk' , Skk' known, one can then easily attain the coefficients matrix of Eg. dTOl) . For Nal, 

2 2 4 

there are totally five linear equations for pp' {Fi = 1) and pp' ' {Fi = 2). By solving them, we obtain 





-2) 












= 2) 



( 

2327^ 
1606 
\39.31-f 



34492- 5511 cos 26*^ + lOOcos40^ 
7289 cos 26*^ - 157.4 cos 46^^ - 23.4 cos 66*, 
h 5823COS26',. - 435.1 cos46'r - 4.5cos66'r 
75.31 cos20r + 35.56cos46'^ + 15.94cos66'ry 



(16) 



where we have defined — Pq{Fi — 1)/(26416— 4549cos20r— 163.4cos46'r+24.5cos66'r). The results are demonstrated 
in FiglH The triangle dependence is caused by the precession of the atoms around magnetic field. As we see, the 
alignment changes sign at the Van Vleck angle Or = 54.7° same as the case for atoms with only fine structures (Paper 
I). As explained in fJJland in Paper I, this is a generic feature of atomic alignment determined by the geometric relation 
of the pumping source and the magnetic field. 

Scattering from such aligned atoms causes polarization in both D lines. Put Eq. lfTO)) into Eqs pTl fc[T^ and combine 
with Eq.(l7]), we obtain the expression for emission coefficients after some tedious calculations. For D2 line, i^; = 1, 



3V3A" 0,, 



vo) {62830 
f(395 



10084 cos 29„ 



- 8.7cos6er + 30.3COS 46^ + 2236.6008 26^ 



412 cos 49,. + 56 cos 69^ + (0.8 cos J 

+ 776) cos 28 + cos20r [(395 - 513.4 cos 29^ + 107.4 cos 49^ + 12.3 cos 69^ - 1.3 cos 89,.) sin^ 
+ cos 4>r [104.3 + 1.95 cos2er - 212.19 cos 49,. - 1.78 cos 69,- + 2.98 cos 89,. - 0.18 cos 109,. 
+ 8.5(cos29r - 1.4 cos 49,. + 0.1 cos 69^ + 48.5) sin 29 sin 29^]} 

(1552.7 + 4473.3 cos 29,. + 60.6 cos 49,. - 17.4 cos 69,. + 1.5 cos 89^) sin^ 9 + cos 20,. [592.5 



3V3A2 0,, 



Sir 

+ cos 29(197.5 - 0.7cos89,. + 6.2cos69,. + 53.7cos49,. - 256.7cos29,. 
+ cos(/)r(17.0cos29r - 23.8cos49r + 1.4 cos 69^) sin 29 sin 29^} 



770.1COS29,. + 161.1COS49,. + 18.5cos69r 



2 cos ^ 



8ir 



-AI,ng°^{i' - va) [cos 9 (790 - 1026.9 cos 29^ + 214.8 cos 49,. + 24.7 cos 69^ + 2.6 cos 89,.) sin 20,. + (34.0 cos 29,. 



* In fact, the absorptions from alkali species are not polarized for the same reason. 
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The numerical coefficients of Eqs. (I1I2I10I [TTIi calculated for NaI according to Egs.©-©. 
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Fig. 3. — Left: The schematic of Na and K hyperfine levels; Middle: Normalized density tensor components Pq/pq{Fi = 1) for ground 
state of Na; Right: for ground state of K. The difference of Na and K alignments is due to different degrees of coherence on the excited 
state (see text). 



F, 



47.6 cos 46^ -1-2.8 cos 69^ + 1649.9) sin <f>r sin 8 sin 29^ 



3^3 A 



(17) 



Al.iigg^iiy - i/o) |l36670 - 28893 cos 28^ + 899cos46l^ - 2cos6er + cos^ 61(20729 cos 28^ - 1118 cos 48^ - 19cos66»r 

;os8e^ + 6057) -I- cos20r (5649.4 - 6278.5 cos 20^ + 656.2 cos 46,. - 27.4 cos 60^ + 0.3 cos 80,.) sin^ + cos 0^ [2522.6 
47.8 cos 20,. - 2521.1 cos 40r + 147.6 cos 60,. - 1.6 cos 80^ -I- 0.1 cos 100,- -I- (5090 - 591 cos 20,. + 6 cos 40,.) sin 20 sin 20,.]} 

^\/^\^ n f o 

AI,ne„i(L' - vo) \ (20729 cos 20,, - 1117.7 cos40,, - 19.3cos60,. - 1.9cos80,. + 6057.1) sin^ + cos 20,, [(-6278.5 cos 20,, 



+ 656.2 cos 40,, - 27.4 cos 60,, + 0.3 cos 80^ + 5649.4) cos^ - 6479.6 cos 20,, + 656.2 cos 40,, 
+ cos0,,(lOO97 - 1183.0 cos 20,, + 12.4 cos 40,, - cos 60,,) sin 20 sin 20,, + 6057.1} 
3V3A 



27.4 cos ( 



. + 0.3 cos 80,, + 5649.4 



— AI,ng?,(,{i/ - i^o) (cos (20189 - 4818 cos 20,, + 215 cos 40,, - 2 cos 60,,) sin^ 0^ sin 20,, 



+ (20193 - 2366 cos 20,, 



-I- 25cos40r - 2cos60rcos80r)sin0,,sin0sin20,,} 



(18) 



For Dl line, 
Fi = l, 



Stt 



AI,ng°^{u - i/q) {36142 - 6350 cos 20^ - 285 cos 40^ + 35 cos 60^ - cos 20(40 - 2 cos 60^ + 76 cos 40^ + 26 cos 20,,) 



- cos20r [(588 - 7O2cos20r + 126 cos 40,, - 11 cos60r) sin 0j - cos 0^ [268 - 1.19 cos 20^ - 267 cos 40^ 
+ cos 60,, + (537.0 - 5.1 cos20r -I- 2.7 cos 40,, - O.3cos60r) sin 20 sin 20^]} 

£1 = _ 3v^A y^j^^^o^f^^ _ I (130.89 cos 20,, - 153.45 cos 40,, + 5.26 cos 60^ + 1.82 cos 80,, - 2050.2) sin^ 
87r L 

-I- cos20r [cos20(-293.92 - 0.20 cos 80,, + 5.66 cos 60,, - 62.56 cos 40^ + 251.03 cos 20,,) + 1053.08 cos 20,, - 187.674 cos 40,, 



8 



Yan & Lazarian 



Q/l of Na D2 emission =90°) U/l of Na D2 emission (^i =0°) Q/l of Na D2 emission (^i =0°) 




Fig. 4. — Contour graphs of polarization signals of Na D2 emission line: left, right: Q/I, middle: U/I. Polarization depends on three 
angles: 9r, 9 and <j)r (Fig[2]l. (f>r is fixed to 7r/2 and 0. At (f>r = 0, U=0. The Stokes parameters Q represents the linear polarization along ei 
minus the linear polarization along 62; U refers to the polarization along (ei + e2)/\/2 minus the linear polarization along (— ei + e2)/\/2 
(see Fig\^ight) . 

+ 16.975 cos 69r - 0.606 cos 861^ - 881.772] - cos 0^ [1073.92 - 10.216 cos 261^ + 5.472 cos 46^ - 0.673 cos 66^] 
sin 26 sin 26^ + 1969.6} 

£2 = _i^^^^i^A7,ne°5(i/ - i/q) [1404.1 cos e(cos2er - 0.1782 cos46lr + 0.0161 cosOe^ - 0.0006 cosSe^ - 0.8373) sin20r 
87r 

+ (20.4 cos 2e^ - 10.9 cos 48^ + 1.3 cos 66,. - 2147.8) sin sin S sin 26^] (19) 

Fi = 2, 



-AI,ng°^{i^ - vq) {64488 - 10148 cos 26>r + 337cos46>,. - 10cos66»^ - cos 26(0.11 cos 86,. + 2.63cos66r - 66.73cos46r + 65.45 cos 26,. 



- 40.301) - cos 20,. ^(-587.8 + 702 cos 26^ - 125.2 cos 46,. + 11 cos 66,.) sin^ ej 

+ cos<^,. [268 - 1.19 cos 26,, - 267.1 cos 46,, + 1.1 cos 66,, - 0.7 cos 86,, + 0.1 cos 106^ + sin 26 sin 26^ 
(536.96 - 5.11cos26,, + 2.74cos46,, - 0.34 cos 66,,)]} , 
ei = -e2(F, = 1), £2 = -£2(F, = 1). (20) 



For unresolved D2 line, the result can be obtained by the summation of the two hyperfine components Fi ~ \ and 
Fi — 2. The angular dependence of the polarization comes from both Pq{Fu, F^), the density matrix of the upper 

level and jQ{i, 51) (see Eq fT2l) . While the former one reflects the direction of incident radiation seen in the theoretical 
frame, the latter is an observational effect which solely determined by line of sight in the theoretical frame. Fig[5] is 
a contour graph showing the dependence of polarization on Or and 9 (see Fig [2]) with fixed (pr = 0°,90°. Along Or 
axis the principle harmonic is of order 2 dependence appears as expected from Eg. (fT7( fT8|) . At 4>r — 90°, Q/I and 
U/I are shifted with respect to each other in Or by a quarter of their period 180°. At 1/)^ = 0°, U=0, and Q/I is 
distorted as the phase dependence of Or and is entangled. Since U=0, the polarization always lies in the single plane 
formed by the radiation source, magnetic field and observer, just as expected from the symmetry of the system (see 
Figl2|). Fig l5|6l are the corresponding plots for degree of polarization p and the positional angle of polarization x (see 
¥\g^ight) calculated according to EqfTil 

Figl?] is the polarization diagram (or Hanle diagram) of Na D2 emission line. Solid lines are the contour of equal 
(pr^ while dash-dot lines are contour of equal Or- For any pair of Or,4'r, the polarization diagram gives the polarization 
Q/I, U/I. The actual diagram is three-dimensional and its shape depend on the perspective we observe (the angle 0, 
see Figl2])- We present here its projection at four directions = 0°, 30°, 60°, 90°. For > 90°, the contour of diagram 
(f>r at 0, is the same as that for 180° — at 180° — 0; and the contour of diagram of Or is the same as that for 180° — Or 
at 180° — 0. This At = 0°, the polarization is symmetric about Or = 90°, i.e. contour of Or coincides with that of 
180° — Or- At = 90°, the polarization is symmetric about Or = 90°. In other cases, the equal-value lines are always 
increasing clock- wise with anti-symmetry of U about 90°. At 6'^ = 0°, 180°, the lines degenerated to a point in every 
diagram (marked by 'x') as expected. At 4>r = 0°, 180°, U = and the direction of polarization traces magnetic field 
in pictorial plane with 90° uncertainty. The lines at (pr = 180° are anti-symmetric in respect to Or = 90° with that of 
(j)r — 0°. These symmetric features are generic and independent of specific species as they are solely determined by 
the scattering geometry. 

The polarizations of the two hyperfine components {Fi — 1,2, see FigEJe/f) of Dl line are reversal to each other 
(see Eq |19l20T) . The polarization of Dl line is thus dependent on the ratio of its line- width and the separation of the 
two hyperfine components ^ Ikm/s. As an illustration, FiglSje/t shows how the polarization of Dl line with a fixed 
radiation and magnetic geometry changes with the line-width. 

Measurements of the polarization degree of both D lines can constrain up to four parameters, from which we may 
extract both magnetic field {Obt4>b) direction and information of radiation source [O^^Wa/W). In future, in the 
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Fig. 5. — Polarizations of Na D2 emission line and their dependence on (9,., 9) at <^r = 90°. Upper panels give the degree of polarization 
of emission line; Lower panels show the positional angle of polarization measured from the plane parallel to magnetic field (see ¥\^^fight). 




Fig. 6. — Degree of polarizations of Na D2 emission line and their dependence on (dr, 6) at 4>r = 0°. The positional angle of polarization 
is zero as the Stokes parameter U=0 (see Fig|4]and the text). 
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Fig. 7. — Polarization diagrams (or Hanle diagrams) of Na D2 emission line observed at different angles d. In the diagram, solid lines 
represent contour of equal (fir, while dash-dot lines refer to contour of equal dr (see Fig[2ll. The numbers in textbox mark the values of 0r, 
while the numbers without boxes give the values of Or- Note that when dr = 0°,&; 180°, the corresponding polarizations degenerate to one 
point in each diagram (marked by o). 
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Fig. 8. — Left: Polarization of Na Dl line at 9r = 90°, = 90°,tt>r = 0°. As we see, the polarization degree and the profile of Dl line 
depends on the line-width. The stokes parameters of the two hyperfine components (Fi = 1,2, see Fig[3{e/i) are reversal to each other. 
Thus if the line- width is much wider than their separation (^ Ikm/s), their polarizations cancel each other out; Middle & right: Contour 
graphs of the emissivity ratio of Na D doublet. 4)r is fixed to 7r/2 and 0. 
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Fig. 9. — The line ratios of Na D2 and DI lines varies with the orientation of magnetic field. Upper: at 0r = 0°; Lower <f)r = 90°. The 
corresponding ratios without magnetic realignment would be the curve of Or = 0° (see Fig[2]l. 

situations when we can resolve the hyperfinc components of D hnes, we can cross-check and make the detection of 
magnetic field even more accurate. 

The intensity of the scattered fight is also modulated by magnetic reafignment. For comparison, we plot the 
D2/D1 line ratios of their intensities (see Yig^niddle & right, FiglS]). The corresponding line ratio without magnetic 
realignment are equal to those values at = 0°, where magnetic field is parallel to the incident radiation. 

4.2. Results for KI 

K I has the same electron configuration and nuclear spin as Na 1 does (see Yig^eft). The only difference between 
them is the coherence between hyperfine sublevels on the excited state 2^3/2- wio — O.5I7, W21 = 1-57, ^^32 = 3.67. 
However, the result for K 1 only differs from that of Na I by < 5%. Here we only provide the density tensors of the 
ground state. 



pI{Fi = 1) 



0.9 cos^" Or + 8.3 cos^ 9r - 44.89 cos^ Br - 135.8 cos'' Or + 535.8 cos^ Or ~ 161.9, 

-2.5 cos^ Or + 15.5 cos*^ Or + 67.0 cos^ Or - 758.4 cos^ Or + 1337.0, | / 

^0.01 cosi2 Or + 0.16 cosio 9r + 6.08 cos* Or + 2.54 cos^ 9,, - 218.68 cos"* 9r + 502.59 cos^ 9r - 143.33, 



/ (-1.0 cos^" 61^ -4.2 cos* ^ 
4,^ ^,-7.4893 COS* 



+ 58.9 cos^ 9r - 57.1 cos'^ Or - 
h 91.6742 cos6 9r - 46.534 cos" 



438.0 cos2 Or + 1026.9) 



(21) 



1.5742 cos-^ 



1.3429 



2.8 cos* 9r - 53.3cos6( 



422.2 COS* 9r + 1320.0 cos^ 9r - 5521.9 



(22) 



There are a number of parameters that determine the polarization: direction of magnetic field (9b,4'b), direction of 
radiation field {9o,4>o) and the percentage of radiation anisotropy Wa/W. For those cases the radiation source is 
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known, (/iq can be easily obtained. If we know the distances of the source and atomic cloud, we can also determine 
9o and Wa/W . Thus we have 2-4 unknown numbers depending on the specific situation. The list of observables we 
have are degree and direction of polarization, line intensity ratio of the doublet. If the hyperfine separation of the 
sodium doublet (~ 20mA) is resolved the four line components are available. In this case one has enough observables 
to constrain the 3D direction of magnetic field and the environment in situ. 

4.3. Comparison with earlier works 

As we mentioned in the Introduction, the optical pumping of atoms have been studied in laboratory in relation 
with early day maser research (see Hawkins 1955, Happer 1972). For instance, magnetic mixing of the occupations 
of Sodium atoms was extensively studied by Hawkins (1955). Their quantitative results, however, is incorrect. This 
is because of the classical approach they adopted for radiation field. As pointed out in Paper I, because conservation 
and transfer of angular momentum is the essence of the problem, it is necessary to quantize the radiation field as the 
atomic states are quantized. Hawkins (1955) claims that the alignment is zero at 6r = 90° {9r is the angle between 
magnetic field and the pumping light). Our results show that the alignment is zero at the Van Vleck angle 54.7°. And 
this is a generic feature for optical pumping of any atoms regardless of their structures (see also Paper I). 

5. TIME-DEPENDENT ALIGNMENT 

In this paper we mostly deal with alignment in equilibrium states. Astrophysical environments may present us with 
the cases when there are not enough scattering events to reach equilibrium. This can happen if either the mean free 
path of the atom is comparable or larger than the system dimension so that atoms are leaving the system before 
acquiring the steady-state alignment, or if the rate of randomizing collisions is comparable to the optical pumping 
rate. In this case, a collisional term should be added to the statistical equilibrium equations (see App lCj) . Below we 
study the alignment of Na I under a given number of scattering events. 

The calculations are straightforward. Initially when no pumping has occurred, the ground state occupation is 
isotropic. There is only Pq in this case. Since the energy splitting between the two hyperfine sublevels is negligible, 
atoms are distributed according to their level statistical weight. 





= 1)\ 


fV3\ 


pUFl 

pI{Fi 
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We simply need to find the explicit expression for scattering matrix and then multiply the initial density matrix by 
it as many times as scattering events. For such purpose, we need to go back to the statistical equilibrium equations 
(|l|2p . The first step, the stimulated emission is described by the second term on the right side of Eq.(IT]). There are two 
routes then for the excited atoms, either they spontaneously decay (first term on the right side of the same equation) 
to the ground level or precess under the magnetic field^ produced by nuclear spin (the 3rd term on the left side). 
The probability of spontaneous emission out of the two routes is thus easily obtained by multiplying the first term on 
the left side of Eq.® by A/ {A + 2TTii/p^^p> ). It turned out that the scattering matrix is the same as the first term 
in Eg. ljlOp apart from the Einstein coefficient B which should be removed as we are concerned about the probability 
rather than rates. The effect of magnetic field is only to remove the magnetic coherence as explained earlier. We 
assume for the incident radiation the same intensity at the resonant frequencies of D lines. Then averaging over both 
D lines, we obtain the scattering matrix 



/1. 09 
0.497 

0.279 cos^ 6 
0.0738 cos^ 

Vo 



0.497 
1.35 

0.0929 0.115 cos^ 
- 0.0246 0.269 cos^ 




0.0340 cos^ - 0.0113 0.052 - 0.156cos^ 9^ 

0.00878 - 0.0264 cos^ 9r 0.121 cos^ Br - 0.0403 
.- 0.0383 0.190 cos^ -(- 0.225 0.105 008^6^ + 0.206 0.0219 - 0.0656 cos^ 6, 
. - 0.0896 0.0879 - 0.0315 cos^ 0.0602 cos^ + 0.699 0.136cos^ - 0.0453 

0.0845 - 0.254cos^ 9^ 0.294cos^ 6^ - 0.0981 0.165 cos^ 8^ + 0.22 



(24) 



For a given number of scattering events, the density matrix of the ground level can be obtained by multiplying 
Eq|23l the scattering matrix Eql24| as many times as the number of scattering events. The results are shown in 
FigllQl We see after ~ 5 scattering events, the system approaches its equilibrium. The scattering time scale is 
BiuJq ~ 5 X 10~ [R{Rq) /r{pc)]'^ . In contrast, the collisional transition rate is ~ 10~^''neCm^ (Happer 1972). In this 
sense, collisions can be neglected if r <C Ipc for a solar type star. In the case where collisions are not negligible, the 
effect of collisions would be to limit the number of scattering events. 

We also obtained the polarization of emission and absorption from atoms scattered one photon. The shape of 
polarization curves for all the lines are very similar to the equilibrium cases. In fact, the curves for positional angle 
of polarization are the same for all the cases. In this sense, direction of polarization is a more robust measure for the 
detection of magnetic field. Only the amplitude of degree of polarization is decreased from the equilibrium case by 
~ %10. 



^ Note, that in this paper we discuss a regime for which the external magnetic field is not strong enough to affect the upper level 
population. 
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Fig. 10. — The density tensors of Na ground state for nonequilibrium cases. The degree of ahgnment pg increases with the number of 
scattering events. Every scattering event is described by an statistically averaged scattering matrix (see text). 

6. EXAMPLE; SYNTHETIC OBSERVATION OF COMET WAKE 

As an illustration, we discuss here a synthetic observation of a comet wake. Though the abundance of sodium in 
comets is very low, its high efficiency of scattering sunlight makes it a good tracer (Thomas 1992). It was suggested by 
Cremonese & FuUe (1999) there are two categories of sodium tails. Apart from the diffuse sodium tail superimposed 
on dust tail, there is also a third narrow tail composed of only neutral sodium and well separated from dust and 
ion tails. This neutral sodium tail is characterized by fast moving atoms from a source inside the nuclear region and 
accelerated by radiation pressure through resonant D line scattering. While for the diffuse tail, sodium are considered 
to be released in situ by dust, it is less clear for the second case. Possibly the fast narrow tail may also originate from 
the rapidly fragmenting dust in the inner coma (Cremonese et al. 2002). 

The gaseous sodium atoms in the comet tail acquires not only momentum, but also angular momentum from the 
solar radiation, i.e. they are aligned. Distant from comets, the Sun can be considered a point source. As shown in 
FigHH the geometry of the scattering is well defined, i.e., the scattering angle Oq is known. The polarization of the 
sodium emission thus provides an exclusive information of the magnetic field in the comet wake. Embedded in Solar 
wind, the magnetic field is turbulent in a comet wake. We take a data cube fFig fT2)) from MHD simulations of a comet 
wake. Depending on its direction, the embedded magnetic field alter the degree of alignment and therefore polarization 
of the light scattered by the aligned atoms. Therefore, fluctuations in the linear polarization are expected from such a 
turbulent field (see Fig ll3ll4| . The calculation is done for the equilibrium case. If otherwise, the result for degree of 
polarization will be slightly different (< %10) depending on the number of scattering events experienced by atoms as 
presented in S}5l The direction of polarization, nevertheless, should be the same as in the equilibrium cases. Except 
from polarization, intensity can also be used as a diagnostic. Note that the result also depends on the inclination with 
the plane of sky a — 90° — Oq . FiglT4] shows that the patterns are completely different at a = 45° from the ones with 
no inclination (Fig ll3p . By comparing observations with it, we can determine whether magnetic field exists and their 
directions. For interplanetary studies, one can investigate not only spatial, but also temporal variations of magnetic 
fields. Since alignment happens at a time scale tji, magnetic field variations on this time scale will be reflected. This 
can allow cost effective way of studying interplanetary magnetic turbulence at different scales. 

7. ALIGNMENT OF H I, PV, & NV 

7.1. Aligned atomic hydrogen 

Hydrogen has a similar structure as sodium does. The nuclear spin of hydrogen is / = 1/2. The total angular 
momentum of the ground state can be Fi = 1/2 ± 1/2 = 0,1 (see FigfTSteft). Only the sublevel of = 1 can 
accommodate alignment. The hyperfine splittings vp^ pi^ oi nPi/2,nP-^/2 are smaller than their natural line width 7, 
V2i{nP^/2) = 0.2297, and vio{nPi/2) = O.2587. Thus coherence on both levels should be taken into account. We 
obtain from Eq. pO|) . 



pUf, = 1) 



(l.SSTcos-*^^ + 3.114cos2 Or - 23.665) 
0.11008- 0.33023 cos2 Or 



(25) 



where g'^ = p[J(F, = 0)/(cos''6'r + 1.730cos2 6'r - 13.652). Insert the ground density matrix Po'^(^/) into Eqs. (|12|13p . 
we obtain their emissivities. For D2 line, 



AqoHi^ - I'd) {cos'^ e(1.092 cos*^ 9^ + 



1.762 cos' 



17. 148 cos'' e,. + 5.481) + cos 20^ sin'' 9(6.131 cos'' 6^ - .3446 cos" 6,.) 





distance *ng the tail distance across the tail distance along the tail ° ° distance across ttie tail 

Fig. 13. — Polarization caused by sodium aligned in the comet wake {a = 0°, Fig llll right). Spatial and temporal fluctuations of 
polarization carry the information on MHD turbulence. Degree and direction of polarization for D2 emission and absorption lines. 
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Fig. 14. — Same as Fig |13l but for the case when comet tail has an inclination (a = 45°, Fig llll right) with the plane of sky. As we see, 
the pattern is completely different for a different inclination. 



2P 



F|=2 
-Fi= 1 



Fi=0 



D2 



Dl 



3S„ 



F,= l 
Fi=0 



Dl 



F,= 5/2 




F,= 3/2 


0.25 


Fi= 1/2 






0.2 


F,= 3/2 




Fi= 1/2 


0.15 




0.1 


cy o 

e 


0.05 









-0.05 


Fi=3/2 


-0.1 







of HI 






of PV . 


\\ 






V\ 

* 

t 

\ 


/' 
// 
'/ 

/- 
l' 
If 
* 
r 




/ 




\ 

\ 


7^ 

♦ 





20 40 60 80 „ 100 120 140 160 180 



Fi= 1/2 



Fig. 15. — Left: The schematic of H I and PV hyperfine levels; Middle: schematic of NV hyperfine levels; Right: The normalized density 
tensor components Pq(Fi = l)/Po(i^i = 1) of ground state of HI (solid lines), PV (dash-dot lines) and NV (dotted lines). As we see, HI is 
less aligned compared to PV and NV. This is due to the high degree overlap of the hyperfine levels on the upper state of H I. 
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For Dl line, 



5388003" 8r + 3.1297 



(27) 



eo 



Ahng"^^{i^ ~ 1^0) (1.749 cos'' 61^ + 3.447 cos2 6'r - 26.292) , ei = 0, £2 = 0. 



(28) 



The polarization of Lyman a lines are given in Fig |16l We see that compared to Na I, K I, Lyman a line is much 
more polarized. In general, the more substate the atoms has, the less their polarizability is. As polarized radiation are 
mostly from those atoms with the largest axial angular momentum, which constitute less percentage in atoms with 
more sublevels. 

The alignment on the ground state also causes the change in optical depth T21 of HI 21cm line, which is a transition 
between the hyperfine sublevels Fi = 0,1 (see FigfTSteft) . 

^21 - |-^™A2e(z^ - lyoHJ^p^oiFi = 0) - IIkJ^pUFi - l)/\/3] 
ovr 
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Fig. 16. — Contour graphs of polarization signals of Lyman a emission: lejt, right: Q/I, middle: U/1. Polarization depends on three 
angles: 9r, and 4'r (Fig[3]|. <f>r is fixed to 7r/2 and 0. At ijir = 0, U=0. 




Fig. 17. — The change of optical depth of HI 21 cm due to alignment. Left: The ratio of the optical depth with alignment taken into 
account r^eal ^^nd the one without alignment counted tq. The occupations on the two hyperfine levels on the ground state = and 
Fi = 1 become very close due to the Lyman cx pumping. That is why taking into account the alignment component pp or not makes a 
big difference. Right: The percentage of effective absorbers of polarized radiation Q of HI 21cm (see Ea l29l l. Note without alignment pp, 
polarization and this quantity would be zeros. If the radiation field has an isotropic component, these results would be reduced by a factor 
Wa/W (see Eo l29l30l and the corresponding text). 

= —AmX^S.ii^ - vo)n 



(29) 



'pI{Fi = 0) - (1 - 1.5 sin^ B)pI{Fi = l)/^/2" 

For comparison, we plot in Fig llTI the ratio of ratio of the optical depth with alignment taken into account Treai and 
the one without alignment counted tq: 

rreai ^ [JoVjj(f; = 0) - ^KJj'p^jFi ^ l)/V3] ^ [- Jj pI{Fi = 1)/n/3] 

[JSpI{Fi^Q)-J^pI{Fi^I)IV^ [J^pI[Fi=0)-J^pI{Fi^I)/V^ 

Since the two hyperfine levels are almost evenly populated in the equilibrium case, the slight change among the 
populations can make a substantial influence on the transmission of 21cm in the medium (see Fig|17p. 

This may be related to the Tiny-Scale Atomic Structures (TSAS) observed in different phases of interstellar gas 
(see Heiles 1997). Besides, as we see from FigfTTl the optical depth can become negative in some cases, indicating the 
amplification of the 21cm radiation, or maser (see Varshalovich 1967). These effects should be included for HI studies, 
e.g. Lyman a clouds (Akerman et al. 2005), local bubbles (Redfield & Linsky 2004), etc. Moreover, polarization can 
occur due to the alignment. The absorption coefficients for the Stokes parameter Q is given by the second_ component 
of Eq.((29|) (see also ¥isi\T7Viaht) . According to Eq. ((30| the real optical depth is approximately oc Pq (x Jq (x Wa/W 
(see Eq l8l9|) . where Wa is the anisotropic part of the dilution factor of radiation field (see Van de Hulst 1950 for its 
expression) . Thus if there is an isotropic component in the radiation field, the results will be reduced by a factor of 
Wa/W. Note if photons are scattered multiple times before reaching the atoms, the anisotropy of the radiation field 
will be diminished and so these effects. Detailed study will be provided elsewhere. 



7.2. Case of PV 
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Fig. 18. — The change of optical depth of NV 70.7mm due to aUgnment. Left: The ratio of the optical depth with alignment taken 
into account r^eal a^nd the one without align men t counted tq; Right: The ratio of the density effective absorbers to the total density for 
polarized radiation Q of NV 70.7mm (see Ea l29l l. Note without alignment p?, polarization and this quantity would be zeros. 

The overlap of hyperfine structure ra upper levels reduces tne alignment. In fact, P V has the same electron 
configuration and nuclear spin (see FigfTSteft) . PV is more aligned as it does not have the overlap on the upper level, 



p'oiFi - 1)' 
pUFi = 1) 



DS^ 6,. + 2 cos2 0,. - 19 
V2(l - 3cos2 0r) , 



(31) 



where Qq — Pq{Fi = 0)/ (1.876 cos"* 9r + 0.289 cos^ 6r — 10.825) . Fis\TE\riaht gives the comparison of the density tensor 
components of H I and PV. 

7.3. Case of NV 

NV is also an alkali atom. The nuclear spin of Nitrogen is / = 1. The total angular momentum of the ground state 
thus can be F = (1 ± 1/2) = 3/2, /2 (see FisHMniddle) . Therefore the ground state has totally (2 x 3/2 + 1) + (2 x 
1/2 + 1) = 6 sublevels which enables alignment. For NV, the hyperfine splitting is much larger than the natural width 
of the excited state. The smallest hyperfine splitting is about 5.67, their influence is thus marginal. We obtain its 
density matrix as follows: 



P°,iFi = 3/2) 
pUF = 3/2) 



_ f 13. 136 cos-' 0^ - 240.06 
^ ^" I 20.90- 62.701 cos^ 6*^ 



(32) 



where = p^{Fi = l/2)/(16.62 cos'' 6*^ - l.Scos^ Or - 167.90). 

The hyperfine transition between the two sublevels on the ground state (A=70.72mm) has been shown a good tracer 
of hot rarefied astrophysical plasmas (Sunyaev & Docenko 2006). Same as HI 21cm, the alignment alters the optical 
depth (according to Eal29p. which one must take into account when analyzing the NV lines. Further more, polarization 
appears as a result of the alignment (see FigfTS)). 

8. MORE COMPLEX ATOMIC SPECIES: N I 

Unlike Alkali atoms, neutral nitrogen is alignable within its fine structure. The ground state of N I is 'iSy2J ^^'^ the 
excited state is 4Pi/2,3/2,5/2- The ground state therefore can have four magnetic sublevels with M = ±1/2, 3/2. Thus 
hyperfine structure is not a prerequisite for alignment. However, the alignment and resulting polarizations will be 
miscalculated if we do not include the hyperfine structure. For resonant lines, hyperfine interactions cause substantial 
precession of electron angular momentum J about total angular momentum F before spontaneous decay. Thus total 
angular momentum F should be considered and the FMp base must be adopted (Walkup, Migdall & Pritchard 1982). 

Nitrogen has a nuclear spin 1=1. Its ground level is thus split into three hyperfine sublevels Fi = Ji — 1, Ji, Ji + 1 — 
1/2,3/2,5/2. The density tensor p{Fi = 5/2) has three components with k — 0,2,4; p{Fi = 3/2) has two components 
with A; = 0, 2; sublevel i^; = 1/2 is not alignable and we only need to consider Pq. Solving Eq. dTOt . we obtain^ 



Po 



■°''(F, =3/2) 



0/ 3 cos^ + 4 cos^ 6*^ - 552 cos'' 6*^- 872 cos^ 6*^ + 12325 ^ 
^ ^0 U cos^ Or - 72.6 cos^ Or - 235.04 cos"* Or + 1887.5 cos^ Or - 600.4^ 



p°'2(F, = 5/2) = 



10 cos^ 9r + 25 cos^ 9r - 765 cos"* 9r - 1017 cos^ dr + 15087 
6 cos^ 9r - 99.2 cos^ 9r - 335.8 cos"* 9r + 4089.6 cos^ 9r - 1322.3 



(33) 



(34) 



° We did the calculation assuming that hyperfine splitting is at least three times the natural line-width and therefore interference term 
is negligible for the excited state. 
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where = pl{Fi ^ l/2)/(0.2 cos^" 6*^ - 2 cos* 6*^ - 8.6 cos^ 6*^ - 217.5 cos'^ 6*^ - 727.6 cos^ 61^ + 8733.1). 

Since N I has a J; = 3/2 > 1, absorption from N I can be polarized unhke Alkah species (see Sjll). For optically thin 
case, the polarization produced by absorption through optical depth t ~ rjad is (see Eg llSp 

Q_^_!n^ 1.5 sin^ 9w% , 1) -A. T (f , , 2)pg {Fi ) 

Ir^ flo I]f^,(-1)^'-^^. [V2T(F,,OK(FO + T(F,,2)p2(^^;)(l_1.5sin2 0)«;2^^J 

where T = < '{} ^ \, Iq is the intensity of background source. We neglect emission here. For a generic case, 



^Fi Fi I 

where the background source (e.g. QSO) is polarized and optical depth is finite, we can obtain in the first order 
approximation 

(/o + go)e-"('+"^/*) + (/o - Qo)e~^^'^''''/""\ 

g= (/o + go)e-"('+"^/"") - (/o - go)e~"(i~''^/*), 

U = Uoe-\V^V„e'\ (36) 

in which Iq, Qq, Uq, Vq are the Stokes parameters of background source. The polarizations of the absorption to Ju — 1/2 
is illustrated in Fig |19l Note that if the incident light is polarized in a different direction with alignment, circular 
polarization can be generated due to dephasing though it is an 2nd order effect. Consider a background source with a 
nonzero Stokes parameter Uq shining upon atoms aligned in Q direction^. The polarization will be precessing around 
the direction of alignment and generate a V component representing a circular polarization 

It ~ Vi Iq iv m h 

where k is the dispersion coefficient, associated with the real part of the refractory index, whose imaginary part 
corresponds to the absorption coefficient t]. -0 is the dispersion profile and ^ is the absorption profile. 
Optical depth also varies with the alignment fFigi20| . The generic expression of the line ratio of a multiplet is given 

by 

ri _ Sf,(-1)^'-^"^ [V2T(F,,0)pjj(fO + T(f;,2)pg(fO(l-1.5sin^g)z^^,„,^ 1 



T2 Ef,(-1)^'-^- {^T{FuQ)pI(Fi)^T{Fi,2)pI{F,){\ - 1.5 sin^ fl^^..2 



(38) 

Similar to the cases without hyperfine structure (Paper I), absorption is determined by only two angles, dr and 9. 
Among them, determines the ground state alignment and dependence occurs from the direction of observation. 
We also see the Van-Vleck effect in FigfTOl Specifically, the polarization is either || or ± to the magnetic field in the 
plane of sky; the switch happens at Van-Vleck angle Qr — 54.7°. 

In Fig ll9|20l we plot together the polarizations and optical depth ratios for N I and S II, result for which is taken 
from Paper I. As we know, N I and S II have the exactly the same term. The difference between them arises from the 
hyperfine structure of N I. In other words, if we do not take into account the hyperfine structure of N I, it would be 
polarized exactly the same way as that of SII. And this can be tested observationally. As we explained above, it is 
usually true that the more complex the structure is, the less polarized the line is. It can also be interpreted by the 
nature of hyperfine interactions. Hyperfine interactions causes the precession of electron angular momentum in the 
field generated by the nuclear spin. Thus similar to the case of magnetic mixing, the polarization is reduced. 

9. DILUTION ALONG LINE OF SIGHT 

For absorption lines, there are atoms far from any pumping source along the line of sight. These atoms are not 
aligned, and we need to take into account the averaging along line of sight. Different components of the atomic 
density matrix are modulated by the pumping differently. For atoms far from a source, their dipole component of 
density matrix p\ is zero. The zero order term /Jq, representing total occupation of a level, however, only changes with 
alignment by < 10%. As a first order approximation, we thus can ignore the variation of pjj due to the pumping and 
adopt a step function for p\ 

- { pl[aligned)M r < Tc ,oq^ 
P^^'^^-y O,forr>rc- ^^^^ 

where is the distance from a pumping source where the collisional transition rate becomes equal to the optical 
pumping rate. In this case, we only need to multiply p\ in Eq. ()29|35|38|) by Na/Ntot^ the ratio of ahgnable column 
density to the total column density along the line of sight. Accordingly the contours fFig fTTll^ do not change apart 
from their amplitudes because the dependencies on Or-^O only appear in the terms containing pg. By overlapping the 
polarization contour maps and the line ratio contour map, we attain the angles Or and 9. Then insert these values into 
either Eg. ([55]) or Eg. ([55]) . we get the ratio Na/Ntot- More precise results can be obtained by making iterations. From 
this ratio, we learn also the conditions in the vicinity of the pumping source. Combining different atomic species, we 

To remind our readers, The Stokes parameters Q represents the Unear polarization along ei minus the linear polarization along 62; U 
refers to the polarization along (ei + e-£)l\/2 minus the linear polarization along {— ei + e'2)/\/2 (see FiglTy-jg/ii) . 
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Fig. 19. — Left-.the schematic of N I hyperfine levels; and the contour graphs of N I (middle) and S II (right) degree of polarization. SII 
does not have nuclea r sp in. The degree of polarization is determined by the dipole component of density matrix <TQ(9r) and the direction 
of observation 6 (Eas l35ll . N I and S II have the same electron configuration. However, N I has a nuclear spin while S II does not. More 
sublevels are allowed in the hyperfine structure of N I, and the alignment is thus reduced and so is the degrees of polarization of N I line 
compared to that of S II line. 



Optical depth ratio (J^=1/2)/(J^=3/2) of N I Optical depth ratio (J^=1/2)/(J^=3/2) of S II 




Fig. 20. — The contour graphs of line ratios of N I and S II. The difference of their ratios is only due to the hyperfine structure of N I. 

can make a tomography of magnetic field as the ratio Na/Ntot varies with each atomic species. As an example, let us 
consider the alignment of H I by an 0-type star (outside H II region) , for which we know the coUisional transition rate 
is Cio/riH = 3.3 x lO^'^'^cm'^s"^. By equating it with the Lyman a pumping rate BJq, we can get Vc — 15pc in cold 
neutral medium where we adopt hh = 30cm'^s^^. If the star is lOOpc away from us, then the dilution factor along line 
of sight would J^^^"^ nndr/ J^^^^"^ nndr ~ 6, where is the size of Stromgren sphere. 

It is also possible that there are multiple independent pumping sources along the line of sight. These situations, 
however, can be easily identified for diffuse interstellar medium. 

10. DISCUSSION 
10.1. Hyperfine splittings 

We have considered alignment of atoms with nuclear spin in this paper. For these atoms, it is the total angular 
momentum, F = J + 1 that should be considered. The prerequisite for alignment in this case is _F > 1/2. There are 
two categories. Atoms like Alkali atoms, Al II, Cu II, etc. would not be alignable without hyperfine structure since 
their electron angular momentum in ground state is J; < 1. Another category of atoms are alignable within fine 
structure, e.g., N I, CI I, II, III, etc. However, calculations of the alignment and polarization would render erroneous 
results if one does not take into account the hyperfine structure of the species. This is because for resonant lines, the 
hyperfine interaction time-scale is shorter than that of resonant scattering. 

The first category of atoms above, i.e. with J; < 1, cannot produce any polarization in absorptions even the ground 
state is aligned within the hyperfine structure. As we explained in SjSl polarization only reduces due to hyperfine 
interaction. For Alkali atoms, the absorption is not polarized in the frame of fine structure. Taking into hyperfine 
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structure does not make a difference in this case to the absorptions. However, the polarization of emission will be 
affected by the alignment. We note, that although the ground level alignment is not a prerequisite for polarization of 
emission for every line, it does affect the polarization of emission. 

We discussed a few examples of elements from our list in Table 1. They were chosen on the basis of astrophysical 
importance as we see it. For instance, sodium D lines are very pronounced ones and easy to measure. It is also 
important for studies of comet wakes, as we discuss in the paper. More calculations should be done in future in 
relation to particular astrophysical objects under study. 

10.2. Polarization of absorption lines 

We studied polarization of absorption resulting from alignment of atoms with nuclear spin and thus with hyperfine 
structure. The degree of polarization is reduced compared to atomic species of the same electron configurations but 
without hyperfine structure. The direction of polarization, however, has the same pattern, namely, either || or _L to 
the magnetic field on the plane of sky (Paper I). And the switch between the two cases happens at the Van-Vleck 
angle Or = 54.7°. In fact, this should be applicable to all absorption lines (including molecular lines) regardless of 
their different structures as long as the following conditions are satisfied. First, pumping light and background light 
are unpolarized; and then it is in the magnetic realignment regime, namely, magnetic precession is faster than the 
photon-excitation rate. 

This fact is very useful practically. It means that even we do not have an exact prediction and precise measurement 
of the degree of polarization of the absorption lines. We can have a 2D mapping of magnetic field on the plane of sky 
(the angle in Figl2lri(7/it) within an accuracy of 90° once we observe their direction of polarizations. In this sense, 
it has some similarity with the Goldreich-Kylafis effect although it deals with radio emission lines. 

For absorption lines, there is inevitably dilution along line of sight, which adds another dependence on the ratio of 
alignable column density and total column density Na/Ntot- For different species, this ratio is different. The ratio 
should be close to 1 for highly ionized species which only exist near radiation sources. The same is true for the 
absorption from metastable state (see Paper I). Combining different species (with different Na/Ntot), it is possible to 
acquire a tomography of the magnetic field in situ. To extend the technique we shall present elsewhere calculations 
for more atoms with metastable states. 

An additional effect that we consider briefly in §6 is the generation of circularly polarized light when linear polarized 
light passes through aligned atoms. This is a new interesting effect the implications of which we intend to explore 
elsewhere. 

10.3. Polarization of emission lines 

In Paper I we dealt with absorption lines of species with fine structure. This paper deals with absorption and 
emission of both the atom species with hyperfine structure and hyperfine plus fine structure. 

The work on emission of atoms with hyperfine structure can be traced back in time. Studies of alignment of neutral 
sodium in laboratory alignment pioneered more that half a century ago by Brossel, Kastler & Winter (1952), Hawkins 
(1955), and Kastler (1956). These experiments revealed that sodium atoms can be efficiently aligned in laboratory 
conditions if atomic beams are subjected to anisotropic resonance radiation flux. However, the calculations that we 
aware of were not satisfactory. For instance, the classical treatment of radiation adopted in Hawkins (1955) does 
not provide the correct measures of the alignment and polarization. Moreover, the calculations were limited by the 
particular geometry of equipment used. On the contrary, our astrophysical applications require us to get predictions 
for arbitrary angles between the line of sight and magnetic field as well as for arbitrary angles between the illumination 
source and magnetic field. 

In the paper we have considered the situation that atoms are subjected to the flow of photon that excite transitions at 
the rate = Bi^Jq which is smaller than the Larmor precession rate ly^, but larger than the rate of disalignment due 
to collisions t~^. For the cold gas with 30 hydrogen atoms per cubic cm, the characteristic range over which the atoms 
can be aligned by an O star is < 15pc for HI due to spin exchange collisions (with rate Ciq/uh = 3.3 x lO^^^cm'^s^^). 

Compared to absorption lines, emission lines are more localized and therefore the dilution along the line of sight can 
be neglected. The disadvantage of emission lines compared to absorption lines is that the direction of the polarization 
of emission lines has a complex dependence on the direction of the magnetic field and the illusion light. Therefore, the 
use of emission lines is more advantageous when combined with other measurements. 

The change of the optical depth is another important consequence of atomic alignment. Such effect can be important 
for HI as was first discussed by Varshalovich (1967). However, the actual calculations that take into account the 
magnetic realignment in ubiquitous astrophysical magnetic fields are done, as far as we know, only in this paper. It 
might happen that the variations of the optical depth caused by alignment can be related to the Tiny-Scale Atomic 
Structures (TSAS) observed in different phases of interstellar gas (see Heiles 1997). Similar effects may be present for 
Lyman alpha clouds and other objects. 

10.4. Implications 

Atomic alignment opens a new channel of information about the physical properties of the various medium, including 
environments of circumstellar regions, AGN and interstellar medium. In particularly, the topology of magnetic fields 
that are so important for these environments can be revealed. What is unique for this new window is the possibility 
of obtaining information about 3D directions of magnetic field. Combining different emission and absorption lines 
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it seems possible to restore the entire strueture of the region, that would include full 3D information, including the 
information about the position of illuminating stars. 

We have done calculations for a number of representative atomic species (see Table 1). These emission and/or 
absorptions from these atomic species, are important lines seen in different astrophysical environments. Indeed, there 
are much more atomic lines that can be studied the same way as we did here. The particular choice of atoms to 
use depends both on the instruments available and the object to be studied. For instance, Sodium D lines have 
been observed in interplanetary medium, including comet tail, and a Jupiter's moon lo. Many lines from alignable 
species have been observed in QSOs and AGNs (see Verner, Barthel & Tytlcr 1994) and they can be used to study 
the magnetic field in situ. Needless to say, this technique can be used to any interstellar medium near an emitting 
source, H II regions, circumstellar discs, supernovae, etc. For intergalactic gas, Lyman a, H I 21 cm, NV and other 
radio lines, one also should be aware of the fact that alignment also changes the optical depth and the emissivity of 
the medium. And with our quantitative predictions of the alignment, one can also attain the information of magnetic 
field in the medium. 

A detailed discussion of alignment of atoms in different conditions corresponding to circumstellar regions (see Diner- 
stein, Sterling & Bowers 2006), AGN (see Kriss 2006), Lyman alpha clouds (see Akerman et al. 2005), interplanetary 
space (see Cremonese et al. 2002), the local bubble (see Redfield & Linsky 2004), etc., will be given elsewhere. Note, 
that for interplanetary studies, one can investigate not only spatial, but also temporal variations of magnetic fields. 
This can allow cost effective way of studying interplanetary magnetic turbulence at different scales. 

Taken together this paper and Paper I provide examples of treating both aligned absorbing and emitting species. For 
the sake of simplicity, we have not considered a few effects that can affect observations. For instance, we did not consider 
effects of the finite telescope resolution for the emission from the regions that have finite curvature of magnetic field 
or random magnetic field component. Such effects arc well known and described in the existing polarimetric literature 
(see Hildebrand 2000). For absorption, however, the effect of telescope finite diagram is negligible if we study the 
absorption from a point source. In addition, we considered emission from optically thin medium, which justifies our 
neglect for the radiative transfer. 

In the present paper we considered the alignment of Nal atoms subjected to a limited number of scattering events. 
However, our approach to describe time-dependent alignment is general and may be easily applied to other species. 
This may be particularly important for studies of transient phenomena using our technique. 

As the resolution and sensitivity of telescopes increases, atomic alignment will be capable to probe the finer structure 
of astrophysical magnetic fields including those in the halo of accretion disks, stellar winds etc. Space-based polarimetry 
should provide a wide variety of species to study magnetic fields with. 

11. SUMMARY 

In this paper we calculated the alignment of various atomic species having hyperfine structure and quantified the 
effect of magnetic fields on alignment. As the result, we obtained linear polarization that is expected for both scattering 
and absorption. We have shown that 

• Atomic alignment of atoms and ions with hyperfine structure of levels happens as the result of interaction of the 
species with anisotropic flow of photons. 

• Atomic alignment affects the polarization state of the scattered photons as well as the polarization state of the 
absorbed photons. The degree of polarization is influenced by mixing caused by Larmor precession of atoms 
in external magnetic field. This allows a new way to study magnetic field direction in diffuse medium using 
polarimetry. 

• The degree of polarization depends on the species under study. Atoms with more levels exhibit, in general, less 
degree of alignment. More importantly, it depends on the angle between the direction of the pumping light and 
the observational direction with respect to the magnetic field embedded in the medium. 

• The direction of polarization depends on the direction of the anisotropic radiation and observation in respect to 
the magnetic field acting upon at atom, if the polarization of scattered light is concerned. 

• The direction of polarization is either parallel or perpendicular to magnetic field, if the polarization of absorbed 
light is studied. The switch between the two options happens at the Van-Vleck angle between the direction of 
magnetic field and the pumping radiation. 

• The intensity ratio of scattered lines or absorption lines are also influenced by magnetic realignment and therefore 
also carry the information about the direction of magnetic field. 

• Absorption and emission of species along the line of sight away from the pumping sources interferes with the 
detected signal, e.g. influence the degree of the measured polarization and the ratio of absorption lines. This 
effect, however, can be accounted for iteratively. 

• If the light incident on the aligned atoms is linearly polarized, as this is a typical case of QSOs, circular polar- 
ization gets present in the transmitted light. 
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• Atomic alignment is an effect that is present for a variety of species and for different terms of the same atom. 
Combining the polarization information as well as using line intensity ratio data allows to improve the precision 
in mapping of magnetic fields and get insight into the environments of the aligned atoms. 

• A steady-state alignment is achieved when after many scattering events. For a limited number of the scattering 
events the alignment depends on this number, i.e. "time-dependent" . While the direction of polarization are 
the same for both cases, the degree of polarization increases with the number of scattering until the steady-state 
alignment is reached. 

• Time variations of magnetic field in interplanetary plasma should result in time variations of degree of polarization 
thus providing a tool for interplanetary turbulence studies. 
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APPENDIX 
RADIATIVE TRANSITIONS 

For spontaneous emission from a hyperfine state JIF'M' to another hyperfine state J' IF, M, the transition proba- 
bility per unit time is 

^4 4 2 2 3 

a = E I < JIFM\V,\J'F'M' > p (Al) 

where V^'" = r • e^'° is the projection of dipole moment along basis vector of radiation, e± — (t^— y)/^/2, Bq — z. 

For hyperfine lines, in the case of weak interaction in which neighboring fine levels do not interact, the electrical 
dipole matrix element for transition from a hyperfine sublevel F' ,M'p of upper level J' to F,Mp of lower level J is 
given by 



Rl,,=< IJFMr\V,\IJ-F-M^ >= (-1)' +'"-■ < JIFnV,\\J-I'F- > 

where Mp is quantum number corresponding to the projection of the total angular momentum F, I corresponds to 
the nuclei spin. The matrix with big " { }" represents 6-j or 9-j symbol, depending on the size of the matrix. When 
more than two angular momentum vectors are coupled, there is more than one way to add them up and form the 
same resultant. 6-j (or 9-j) symbol appears in this case as a recoupling coefficient describing transformations between 
different coupling schemes. 

IRREDUCIBLE DENSITY MATRIX 

We adopt the irreducible tensorial formalism for performing the calculations (see also Paper I) . The relation between 
irreducible tensor and the standard density matrix of atoms is 

p^{F,F') = J2 i-ir-'\2K+l)'/' ( 2 M' ) < F^'I\P\P'^^I' > ■ (Bl) 

MM' ^ ^ / 

For photons, their generic expression of irreducible spherical tensor is: 



qq' 

EFFECT OF COLLISIONS 



Collisions can cause transitions among the hyperfine sublevels and reduce the ground state alignment. In the regime 
where collisions are not negligible, the equilibrium equation for the ground state Eqs.([2]) should be modified to include 
the collisions. For the ground hyperfine level, 

Pt{F?)= E pM]C{FI ^ F^)pl{Fl) + E Pk[Ju]A{Ju^ Ji)p\{F^,Fi) 



E 6kk'CiF^ ^ FD + D'^ + J2 (Skk'BiJ^ + Skk'BiJ^) 



PoiFn, (CI) 
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For other hyperfine levels on the ground state, 

f;\e>e' -f/Ib'>b J.,f„,_f; 



E Skk'CiFi^Fl)+D''+ E 



poiPi)+ E [JiWciFi ^ FOpUm) 

FI\e'<e 



where C is the coUisional transition rate among the hyperfine levels on the ground state, D'^ is the depolarizing rate 
due to elastic collisions (see Landi Degl'Innocenti & Landolfi 2004). E and E' are the energies of hyperfine levels Fi 
and F/ respectively. 

FROM OBSERVATIONAL FRAME TO MAGNETIC FRAME 

In real observations, the line of sight is fixed, and the direction of the magnetic field is unknown. Thus a transforma- 
tion is needed from the observational frame to the theoretical frame where the magnetic field is the quantization axis. 
This can be done by two Euler rotations as illustrated in Fig[21 In the original observational coordinate (xyz) system, 
the direction of radiation is defined as the z axis, and the direction of magnetic field is characterized by polar angles 0b 
and ^B- First we rotate the whole system by an angle (j)B about the z-axis, so as to form the second coordinate system 
x'y'z' . The second rotation is from the z{z') axis to the z" axis by an angle 9b about the y'-axis. Mathematically, the 
two rotations can be fulfilled by multiplying rotation matrices. 



cos 9b — sint 

1 
sin6'B cos 9, 



cos (f>B sin^sO 
— sin (pB cos 0s 
1 



cos 9b cos (j)B COS 9b sin 0b — sin 9b 





— sin (pB 



COS (pB 



sin 9 B cos (pB sin 9 b sin (f>B cos9b 



(Dl) 
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